Non-perturbative Dynamical Decoupling Control: A Spin Chain Model 
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This paper considers a spin chain model by numerically solving the exact model to explore the 
non-perturbative dynamical decoupling regime, where an important issue arises recently []J. Our 
study has revealed a few universal features of non-perturbative dynamical control irrespective of the 
types of environments and system-environment couplings. We have shown that, for the spin chain 
model, there is a threshold and a large pulse parameter region where the effective dynamical control 
can be implemented, in contrast to the perturbative decoupling schemes where the permissible 
parameters are represented by a point or converge to a very small subset in the large parameter region 
admitted by our non-perturbative approach. An important implication of the non-perturbative 
approach is its flexibility in implementing the dynamical control scheme in a experimental setup. 
Our findings have exhibited several interesting features of the non-perturbative regimes such as 
the chain-size independence, pulse strength upper-bound, noncontinuous valid parameter regions, 
etc. Furthermore, we find that our non-perturbative scheme is robust against randomness in model 
fabrication and time-dependent random noise. 
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I. INTRODUCTION 



Central to quantum science and technology is to com- 
bat decoherence caused by inevitable external noises or 
quantum operation inaccuracies [2|. Strategies for con- 
trolling a quantum state include the quantum error cor- 
rection codes [3H5} and dynamical coupling control such 
as fast-strong pulses ("bang-bang") control |6l-[l2|. One 
common feature of these dynamical control approaches 
is that almost all the theories are defined in a perturba- 
tive manner. Despite the progress made in theoretical 
studies and experimental realizations, the perturbation 
theories still lack a satisfactory understanding of their 
validity and applicability, especially the divergence is- 
sue due to the onset of a possible phase transition. Re- 
cently, a careful examination of the dynamical decou- 
pling or bang-bang control has led to new breakthroughs 
achieved through introducing a non-perturbative dynam- 
ical decoupling scheme based on exact stochastic master 
equations in bosonic baths [l[. By using this scheme, we 
showed that there is a threshold and a large pulse param- 
eter region where the effective dynamical control can be 
implemented, in contrast to the perturbative decoupling 
schemes where the permissible parameters converge to a 
point or form a very small subset in the large parameter 
region admitted by our non-perturbative approach. Our 
non-perturbative approach has shown impressive flexi- 
bility in implementing the dynamical control scheme in 
experimentation. 

This paper uses a spin chain model to explore the non- 
perturbative dynamical decoupling regime. Our results 



demonstrate, through new type of environment and new 
numerical methods, several universal features of the non- 
perturbative dynamical decoupling emerge. More specif- 
ically, we show that there is a threshold and a large pulse 
parameter region where the effective dynamical control 
can be implemented. 

We model the first spin in an XY-type model as our 
system of interest and the rest as our environment. The 
bath spectral density is chosen as a complicated function 
indicating that the ratios between two arbitrary levels 
are normally not rational numbers. We have derived an 
exact "master equation" for our spin system. It should 
be noted that for the spin environment, the environmen- 
tal correlation function is very different from that for a 
bosonic model. However, our results obtained from this 
particular XY-type model are expected to be applicable 
to other spin-environment models. Our numerical calcu- 
lations show that effective control can be made in many 
different ways either through a multiple pulse sequence, 
as in the case of bang-bang control, or just a single pulse 
applied within the bath memory time. The effectiveness 
of the dynamical control is not sensitive to the size of 
spin chain, an observation that showcases the universal- 
ity of our findings. For the spin chain model, we have 
shown that there is an upper-bound for pulse strength in 
non-perturbative regimes when the pulse duration is suf- 
ficiently large, which is rather counter-intuitive from the 
point of view of the strong fast pulse control (bang-bang). 
We also consider the time-independent randomness for 
coupling constants and site energies of the XY-type spin 
chain. We show that our results and conclusion are ro- 
bust to the random couplings and site energies even in 



the strong-coupling limit. Finally, we consider the time- 
dependent random noise, which may simulate the noisy 
effects due to additional environmental variables. We 
show that the quality of dynamical control in the non- 
perturbative regime is also robust to the influence of the 
time-dependent noise. 



II. THE MODEL 

The spin chain model considered here is represented 
by the following Hamiltonian, 



H = H S + H B + H sb 



(1) 



where Hs and Hb is the system's and environment's 
Hamiltonian. Hsb is the interaction between the system 
and environment. We consider an XY-type spin-chain 
model, where the first spin of the chain is our system 
of interest and the environment is composed of all the 
other N — 1 spins. The system Hamiltonian and control 
are Hs = c(t)Z\, where c(t) is a time-dependent control 
function. The control is only applied to the system, and 
the c(t) is chosen as a sequence of periodic rectangular 
pulses: 



it) 



^ nr < t < tit + A, n is integer 
0, otherwise 



(2) 



where ^ is the pulse strength, 1/r is the pulse frequency 
and A is the width of the pulse in one period. The 
three parameters span the parameter space of the non- 
perturbative dynamical decoupling scheme as in PJ ■ The 
limiting case with 1/r, A — > and \I/ — > oo corresponds 
to the idealized Bang-Bang control. 

The XY-type spin-chain Hamiltonian is modeled as 



Hsb = -JiMXiX 2 + Y 1 Y 2 ), 



(3) 



N-l 



Hb — — / J Ji,i+l{XjXi + i + YjYi+x) + -ffsite (4) 

1=2 

where Jj^+i is the exchange interaction between sites i 
and i + 1. The system and environment coupling constant 
is Ji.2- Xi,Yi denote the Pauli operators acting on spin 
i, and N is the total number of sites. We consider a 
natural configuration of the spin chain with open ends. 
The ^-component of the total spin is a conserved quantity 
implying the conservation of the total excitations in the 
system. For simplicity, we consider that the system is 
in the "one-magnon" state, in which the total number 
of up spins is one and J = Jj,i+i = 1.0. -ffsite is the 
site energy and will be specified in our later discussion. 
Our numerical calculations show that, within the "one- 
magnon" subspace, the sign of J (plus or minus) does 
not affect the results. 

The system dynamics alternates between control and 
free evolution. We start with a control pulse with the 
width A, followed by r — A free evolution. The evolution 



operator of the whole system for the first time interval or 
period will be U (t-A)U(A), where 17(A) = cxp{-iAH] 
and Uo(t — A) = cxp{—i(j — A)H ]. Here H denotes the 
free evolution without control. We then repeat the same 
operations from time r to 2r, • • • . The bang-bang control 
theory is the idealized limit, A, r — > and "J — > oo, which 
has been shown to be able to eliminate the interactions 
between the system and bath. 

We can diagonalize the Hamiltonian H and Hq such 
that H d = W^HW, H od = V^HqV. The evolution oper- 
ators can therefore be expressed by 

[7(A) = W cxp[-iAH d ]W f (5) 

Uq(t-A) = V exp[-i(r - &)H M ]V* (6) 

Suppose that our system, the first spin, is initially in the 
state |0(O)) = |1), while the other spins are in the spin 
down state |0). The spin chain will be in a product state 
|$(0)) = |1> <g> |00 - - - 0> = |100- ■ -0>. 

After m-time control, the fidelity at the time mr, mea- 
suring the survival probability of the initial state |^(0)), 
can be defined as F = ^/(</>(0)| p(t) |</>(0)), where p{t) is 
the reduced density matrix of the state at the first sites 
at t — mr. The N eigenvectors and eigenvalues of the 
Hamiltonian iJ, Hq can be obtained by numerical diago- 
nalization. 

In the next section, we will analyzes the quality of 
dynamical control in terms of the fidelity defined above. 



III. RESULTS AND DISCUSSIONS 




FIG. 1: (Color on line) The fidelity after 128-time control 
pulse (t = 128t), the initial state is |1) at the first spin. The 
bottom triangular area has no physical meaning since r is 
always larger than A by definition. W = 8.0, ./V = 130. 

Suppose that our system, the first spin, is initially in 
the state |1). With time evolution, the state will be dis- 
tributed to the other spins due to the interaction among 



the spin chain. The periodical rectangle pulses will sup- 
press the disturbing process, a well-known fact in the case 
of bang- bang control. Fig. Q] plots the contour fidelity as 
a function of the pulse period r and pulse duration A at 
time t — 128r (the A — r phase diagram). It clearly shows 
that, similar to the dissipative bosonic bath [l|, there is a 
large parameter region where the dynamical decoupling 
(or noise suppression) works equally well as shown in 
the red zone in Fig. [TJ The left-bottom corner(A and 
t —> 0) represent the idealized pulses dynamical decou- 
pling, which has the same fidelity as the regimes with 
parameters t and A in the red zone. Interestingly, there 
is an upper-bound for the parameter A « 1.3 in this spin 
model. The onset of the failure of an effective control is 
the fundamental restriction of the bath memory time to 
be discussed later. 



affect our results and conclusions. 
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FIG. 2: (Color on line) The fidelity as a function of number 
of sites N under free and control Hamiltonian, where $ = 
8.0, A = 1.2, r = 1.3 and t = 128r. 

The analysis presented above is for the case when the 
number of sites, N—130. Since our environment con- 
tains a finite number of spins, it is important to analyze 
the length dependence of our control scheme, in order to 
eliminate the finite size effects. In Fig. [51 we plot the fi- 
delity of the system state at the time t — 128r for varied 
chain lengths or the numbers of sites. We compare the 
fidelity of free evolution with the case that the dynami- 
cal control is applied in this figure. The horizontal axes 
represents the number of sites N and the vertical is the 
fidelity. There are several interesting new features in this 
figure. First, for the free evolution, while fidelity is not 
good as expected, clearly the fidelity shows significant 
oscillation for odd-even numbers of sites, whose magni- 
tude decreases with N. This is a typical finite-size effect. 
However, the quality of the non-perturbative dynamical 
decoupling, in the non-perturbative control pulse region 
with A = 1.2, r = 1.3, does not depend on the chain 
length N, as shown clearly in the figure. The fidelity is 
close to 0.98 for all the total numbers of spins N. This 
implies that the finite size effect of the chain, dose not 



FIG. 3: (Color on line) The fidelity as a function of control 
times t. * = 8.0, A = 1.2, r = 1.3, N = 130. We explain these 
curves in the text. 

The above analysis is for the case with a fixed time 
t = 128t. Now we discuss the time evolution of the 
fidelity. Fig. [3] plots the fidelity versus time t. We com- 
pare the five cases: free evolution, and controlled evo- 
lution under the one of the following four conditions: 
(1) constant coupling J, (2) band broadening, (3) ran- 
domness coupling J + 7rand(i) and (4) stochastic cou- 
pling J + ?7rand(i). Here the band broadening refers to 
the uncorrelated random potentials, which is defined by 
-Hsitc = eX)i ranc K*) ( e — 0-5 hi the figure) and gener- 
ated in the process of physical fabrication. 7rand(z) is a 
random function of sites i, where 7 (taken as 0.5 in the 
figure) is the magnitude of the random function, rand(i) 
denotes a stochastic number uniformly distributed in the 
interval (—1,1). It is noticeable that the dynamical con- 
trol quality, characterized by the fidelity, does not change 
significantly with these two types of random distribu- 
tions, even in the long-time limit, which indicates that 
our non-perturbative dynamical decoupling control is ro- 
bust against the defect or frustration in experimental fab- 
rication of the spin systems. Once the spin-chain is gen- 
erated, it does not change with time. However, it is im- 
portant to consider a time-dependent random noise due 
to an external agent (or environment). Equally interest- 
ing, our numerical scheme allows us to model the noise as 
the perturbation by r]rand(i) (rj = 0.1 in Fig. [3]), which 
changes randomly before each pulse period, i.e. the ran- 
dom number 7/rand(i) is constant within one period r 
but varies for different periods. Surprisingly, the control 
quality is shown to be robust compared to the cases with- 
out the time-dependent noises, as shown by the inset of 
Fig. [3J Again, our analysis shows that another univer- 
sal features of our finding based on the non-perturbative 
dynamical decoupling scheme. 

In Fig. HJ we plot the contour fidelity as a function of 
pulse strength \1/ and r/A for given A = 0.1,0.5, 1.2, re- 
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is only a lower-bound. In other words, only bang-bang re- 
quires strong pulse strength, while non-perturbative con- 
trol does not. 



IV. ENVIRONMENT CORRELATION BY 

FESHBACH PROJECTION-OPERATOR 

PARTITIONING TECHNIQUE 

Using PQ partitioning technique [1] , an n-dimensional 
wave function ip t can be divided into two parts: an inter- 
ested one-dimensional vector P(t) and the rest (n — 1)- 
dimensional vector Q(t). i/j t ,H e ff can be simply written 
as 



A 



P 

Q 



H eff 



h R 

R T D 



(7) 



where the lxl matrix h and (n — 1) x (n — 1) matrix 
D are the self-Hamiltonian living in the subspaces of P 
and Q. For our model, h is the control function c(t) and 
R = [ J ... ] is 1 x (n — l)-dimensional. 



D = 



- J 
J 


. 0" 
. 




. o 


. J 
. J 0. 



(8) 



10 12 14 



In the selected one dimensional subspace, P satisfies 




FIG. 4: (Color on line) r/A — $ phase diagram of the fidelity 
at a control times m=128. N = 130. (a) A = 0.1, (b) A = 0.5, 
(c) A = 1.2. 



spectively. There is always a region where the dynamical 
control works with high fidelity, as shown in the red zone 
where F > 0.95. Surprisingly, there are a lower-bound 
and an upper-bound for the strength in non-perturbative 
regimes. For instance, when A = 0.5, this region is 
6.0 < * < 9.0 when 1.0 < r/A < 1.6. When the param- 
eter is close to the bang-bang regime, e.g., A = 0.1, there 



P = hP-i / dsg(t, s)P(s) 



(9) 



where in our case g(t,s) — g(t — s) = Re~ tD ( t ~ s ^R T 
RLe~ ld ^~ s ^ L^ R T , d is the diagonalized matrix, d 
L^DL. Specifically, 



g(t-s) = J 2 J2\ L i 



fc , a e-«K*-> 



(10) 



where L±k are matrix elements of the first row and the 
kth column of the matrix L and E^ is the A:th eigenvalues 
of D. 

In general, Eq. Q is a one-dimensional exact equa- 
tion [13| for arbitrary baths and has to be numerically 
solved. However, g(t,s), is often reduced to g(t — s), is 
a bath correlation function and often has common be- 
haviors that hardly depend on the details of baths. Nor- 
mally, the real part of the function g(t — s) starts from 
one and decays (rapidly) to zero with time t — s, unless 
baths have limited numbers of frequencies, for instance, 
the spectral distribution of a few Harmonic oscillators. 
The imaginary part corresponds to an oscillating factor 
and is not responsible for the leakage from the P space, 
resulting in the renormalized energies. Fig. [S] shows the 
behaviors of the real part of the environment correlation 
function of this specific model. The spectrum structure 
of D in our model is sufficiently complicated such as the 
ratios of any two energy levels are not rational numbers. 




free evolution 
randomness 1 
randomness 2 
randomness 3 
randomness 4 



FIG. 5: (Color on line) The real part of bath correlation g 
versus t-s, where N = 130 and J — 1. The different dashed 
curves denote different random distributions. 



We see that the correlation function (memory function) 
vanishes at the lifetime « 1.7, with a few small damp- 
ing oscillations due to the finite number of spins. Fig. [5] 
also shows that the effect of randomness on the correla- 
tion function (7rand(i)) of J (See the four dashed lines). 
The randomness only affects the correlation function af- 
ter lifetime and different kinds of randomness yield the 
out-of phase curves, meaning that the randomness may 
eliminate the possibility of the revival of the memory 
after the lifetime. Strong and fast pulses suppress the 
decay rate and effectively prolong the lifetime of the cor- 
relation function. This is one important mechanism un- 
derlying the dynamical decoupling control (see, e. g., 
|12j | and references therein). On the contrary, our valid 
parameter region may allow very wide and slow pulses, 
e.g. t « 1.3 (almost the lifetime), which is obviously 



beyond the known perturbative suppression mechanism. 
Our calculations confirm further that the valid parame- 
ter r should be chosen within the lifetime. That is, the 
dynamical control is effective when the environment is 
still in non-Markovian regimes. 



V. CONCLUSIONS 

In conclusion, we have investigated a non-perturbativc 
approach to dynamical control theory through an XY- 
type spin model. We have shown that dynamical control 
with high fidelity can be achieved in a large control pa- 
rameter region. There are many interesting features aris- 
ing from this spin-chain model including the chain-size 
independence, the pulse strength upper-bound, noncon- 
tinuous valid parameter regions, etc. 

Our initiative has opened the new avenue for exploring 
the dynamical decoupling scheme in a more realistic situ- 
ation allowing non-idealized pushes. Our results suggest 
that the non-perturbative approach may be applicable 
to various physical systems such as the liquid state NMR 
with a small number of spins |14l - ll6l |. solid-state NMR 
]YA LL8J , optical lattices [lj| , or quantum dots [20] . 
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